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The electrical dc-conductivity of disordered, non-interacting electrons is calculated in the asymp- 
totic limit of high lattice dimensions d oo. To go beyond the lowest-order contribution in the ex- 
pansion parameter 1/d of the single bubble diagram, vertex corrections are calculated from an asymp- 
totic expression for the two-particle vertex. A mean-field approximation for the dc-conductivity 
containing the leading high-dimensional vertex corrections is proposed which is free of spurious 
non-analyticities, i.e. the conductivity is non-negative and shows no unphysical behavior in d > 3. 



Electronic problems with interactions or disorder can almost never be solved exactly, except for special limits. The 
resolvent functions of an exact solution fulfill conservation laws and have the correct analytic (Herglotz) properties, i.e. 
do not contain spurious poles. This is generally not true for approximate, e.g. perturbative, solutions. Only global, 
self-consistent approximations have a chance to be free of unphysical behavior and to yield the desired analyticity of 
a solution for all input parameters. 

For non-interacting tight-binding electrons in a random potential the fujs± self-consistent solution with the correct 
analytic properties was the "coherent potential apprjoximation" (CPA).El'El The Herglotz analyticity of the CPA- 
equations was proved explicitly by Miiller-Hartmann.H Only later the CPA was found to be the exact solution of 
the Anderson disorder model in two particular limits. First, the CPA with a semicircular disorder distribution was 
shown to correspond to the exact solution of an n-orbital model in theJ-imit n — oo.u Then, after the limit of high 
spatial dimensions d for fermionic lattice models had been introduced,El it was found that the CPA represents the 
exact solution of the Anderson disorder model in d = cx) for arbitrary disorder distributions.EfLl Since then the limit 
d oo has served as a useful tool for deriving self-consistent, fully dynamical approximation schemes for interacting 
lattice electron systems,u referred to as Dynamical Mean-Field Theory (DMFT). 

By calculating a physical quantity in d = oo one obtains a particular mean-field value. The situation becomes 
subtle if the value obtained in this limit is zero. This is, for example, the case for non-local quantities such as the 
off-diagonal propagator Gij, i j. They depend on the distance between two or more different lattice sites (i.e., 
their Fourier transform is wave-vector dependent) and are thus necessarily proportional to some power of e.g., 
Gij ~ 0{1/Vd), for nearest-neighbor sites However, that does not imply that these quantities can be neglected 
in d = oo. They may contribute, since they appear in lattice sums where the summation over the d — > oo many sites 
compensates their 1/d-smallness. To include these quantities properly, one has to calculate their asymptotic behavior 
in the limit d — > oo, thus going beyond the strict d = oo limit. The necessity to go beyond d = oo is also evident from 
the fact that higher-order non-local Green functions are related to lower-order local Green functions by fupctional 
derivatives vLa-generalized Ward identities, reflecting conservation laws. It was recently shown by one of uaj and by 
Hettler et al.ll3 that, when only the leading asymptotic contributions to one- and two-particle Green functions are 
taken into account, the Ward identities are not fulfilled in d = oo. One has to go, beyond the leading order at the 
lower particle level, or introduce anomalous functions, to restore conservation laws.l3 This shows that the definition of 
two- and higher-order Green functions is ambiguous in the limit d = oo. 

A particularly important example of a quantity whose mean-field value in d = oo vanishes, is the electrical 
conductivity a. It is defined from a Kubo formula with the current-current correlation function. In the limit d — *■ oo 
the optical dc-conductivity is given by a single bubble diagram, with a ~ 0( 1/d) .til In a formal 1/d-expansion this 
result would be the first non-vanishing contcibution to a. It is useful to consider this non-vanishing asymptotic result 
as a "mean- field value" of the conductivityli3. Likewise one may define a mean- field value of any physical quantity 
through its leading non- vanishing asymptotic result in the limit d — > oo. 

A mean-field result for the dc-conductivity a ~ 0(l/d) defined in this way does not contain vertex corrections. 
Hence it does not include the physics of back-scatterings. However, in random systems vertex corrections are known 
to be extremely important since they are responsible for Anderson localization at zero temperature in sufficiently low 
dimensions (d — 1, 2) or for sufficiently strong disorder in 3 < d < oo. At least from a diagrammatic point of view it 
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is not yet fully understood how the conductivity a approaches zero at the localization transition.EjO Clearly one has 
to go beyond the mean-field single-site diagrams to incorporate localization effects. 

Most recently Jarrell and Krishnamurthytj introduced systematic non-local corrections to the CPA on the one- 
particle level using the Dynamical Cluster Approximation to obtain results compatible with Herglotz analyticity, 
i.e., a non-negative density of states. Here we choose another route to go beyond the mean- field limit and calculate 
non-local corrections to the CPA two-particle irreducible vertex function. The aim of our paper is to employ the limit 
of high lattice dimensions tojmprove upon the mean-field conductivity a ~ 0(1/ d) by including vertex corrections. 
We follow the proposal of refJj where the high-dimensional asymptotics of the full vertex function with leading vertex 
corrections to the electrical conductivity was derived. Our main result in this paper is a mean-field expression for the 
electrical conductivity which includes leading asymptotic vertex corrections in d ^ oo while remaining non-negative 
in d > 3. _ 

The vertex function from ref.u contains the leading 1/d-asymptotics of all two-particle quantities. Employing the 
Kubo formula for the electrical conductivity with the two-particle vertex we may derive the leading asymptotics for 
the conductivity and its vertex corrections. However, approximations of the full vertex function in the Kubo formula 
can, in principle, lead to unphysical results. Indeed, the vertex corrections to the single-bubble term may have a 
negative sign and hence positivity of the conductivity cannot be warranted. Clearly, a meaningful and consistent 
approximation for the conductivity must never become negative. To obtain such an approximation we represent the 
full vertex function by means of an irreducible vertex and a Bethe-Salpeter equation in the electron-hole channel. To 
produce the leading contribution to the vertex corrections in d ^ oo, the irreducible function must be evaluated with 
its leading non-local contribution. By simplifying the Bethe-Salpeter equation in high spatial dimensions we then 
obtain a closed, mean-field expression for the conductivity with vertex corrections. 



In the following we consider the Anderson disorder Hamiltonian 
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to describe the effects of randomness. Here t* is the hopping matrix element between nearest neighbors, scaled in such 
a way as to produce a meaningful limit d oo, with Z as the coordination number of the latticed The local, static 
potential Vi is a random variable with site-independent distribution function. The conductivity of a quenched random 
system without interparticle interactions is described by averaged one- and two-particle Green functions (resolvents) 
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respectively. It is our first 



goal to determine these functions in the asymptotic limit d oo. 

It is straightforward to derive the d ^ oo limit of the self-energy which carries the information about how the 
randomness infiuences the motion of a single electron. The self-energy becomes local and can be obtained from the 
single-site equation 
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(2) 
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This is precisely the well-known CPA equatioifl for the self-energy. Here the local (diagonal) one-particle propagator 
is denoted by G{z) ~ N^^ J^k G(k, z) = J dp{e)[z — E(z) — e]^^ where p is the density of states. It is less evident how 
to derive expressions for averaged two-particle functions which are consistent with the local self-energy, since we have 
to work explicitly with non-local quantities. At the two-particle level one has to keep two separate lattice points to 
derive the leading asymptotics for large da It is more convenient and practical to work with an averaged cumulant, 
or better with a vertex P defined in momentum space as 



P(ki,Zi,k2,Z2;q) =G-l(ki,Zi)G-l(k2,Z2) [G(2)(ki,Zi,k2,Z2;q)-5(q)G(ki,Zi)G(k2,Z2) 

xG-i(ki +q,zi)G-i(k2 +q,z2) ■ 



(3) 



It was shown in ref.B that the asymptotic d ^ oo solution for the vertex P can be represented as a sum of three 
contributions, solutions of Bethe-Salpeter equations in three inequivalent two-particle irreducibility channels. The 
integral kernels of these equations are always the local two-particle irreducible vertex (being the same in all channels) 
accompanied by a non-local two-particle bubble. The local two-particle irreducible vertex in d = oo reads 



A(Z1,Z2) = 



1 



S^zi) ^ 

SG{z2) G(zi)G(z2) 



1 - 



1 



1 



1 + (E(z) - V^) G(zi) 1 + (E(^) - V^) G(z2) 



(4) 



2 



and the two-particle bubbles containing the entire momentum dependence are given by x^(q;zi,Z2) = 
N~'^ G(k, 2i)G( q±k, Z2). In the limit d oolhe vertex functions from the three two-particle channels (electron- 
hole, electron-electron, and vertical) take the formEj 

r^^(q;^i,^2) 
r"(q;^i,^2) 

The local part of the vertex functions in (|^) is always the same, i.e. is given by 7(^1,22) = 
A(zi, Z2)/ [1 — A(zi, Z2)G(zi)G'(z2)]. The full vertex is a sum of the above three contributions, where the trans- 
ferred momentum q has a different meaning in each channel. This is due to the fact that the irreducibility channels 
are topologically inequivalent and differ in the momentum that is conserved during multiple scatterings. If the in- 
coming particle and hole carry momenta ki and k2 then the conserved motaentum is k2 — ki, ki + k2 + q, and 
q for the electron-hole, electron-electron, and vertical channels, respectively.!!!! The momentum q is the momentum 
transferred during the scattering on impurities, i.e., the outgoing particle and hole carry momenta ki -|-q and k2 -l-q, 
respectively. In order to avoid multiple summation on the same site we must subtract the local vertex from the sum 
of the channel-dependent vertex functions. twice. We then obtain an explicit representation for the two-particle vertex 
in high dimensions in the notation of rcf.!3 

r(ki, zi, k2, Z2; q) = r^^ka - ki; zi, zj) + ^"(ki + k2 + q; 21,^2) + r''(q; 21,22) - 27(21, 22) . (6) 

We note that the CPA vertex function derived in ref§ is given by only the first term in the above equation, i.e. . 
Hence it does not contain the transferred momentum q needed to incorporate vertex corrections to the conductivity. 

The density of the static (dc) electrical conductivity at zero temperature is defined by a Kubo formula with the 
full vertex as {h— 1) 

I^e = |^^5]z;o(k),;;3(k')^(-aT)G,(k)G,(k) [5(k - k') + r,.(k,k;k' - k)G,(k')G,(k')] (7) 

k,k' err 

where cr,T = ±1, r^^(k, k';q) = r{k,EF +iaO+,k',EF -t- irO+; q), G„(k) = G(k, + iaO+), Vc,{k) = 
m~^de(k)/dka, e(k) is the dispersion relation and m the mass of the electron. 

Eq. (0) with the vertex functions (|^) contains nontrivial corrections to the one-electron conductivity (the single- 
bubble diagram). However, (|^) is not appropriate for approximate evaluations in finite dimensions. The vertex 
corrections are merely added to the one-particle conductivity, such that negative contributions may reverse the overall 
sign, thus leading to unphysical behavior. To avoid such a situation we represent the conductivity in a different way. 
We use a Bethe-Salpeter equation in the electron-hole channel expressing the full vertex F via an irreducible one, A'^''. 
The irreducible vertex A^'^ together with x+ determine the integral kernel of the Bethe-Salpeter equation explicitly 
and define a matrix multiplication scheme in momentum space.!il! The integral kernel and the multiplication rule for 
the electron-hole channel are given by 

[Af^G^G,] (k,k';q) = Af,(k,k';q)G,(k-hq)G,(k'+q) , (8a) 

[X . Y] (k, k'; q) = ^ E ^(J^' 1^'' q')>^(k + q', k' + q'; q - q') . (8b) 

q' 

A solution of the Bethe-Salpeter equation for the two-particle vertex A'^'' can formally be written as F = 

{1 - [A'^'^GG] •}" A'='' where the bullet indicates that, upon expansion of {...} ^, the two-particle functions [A'^'^GG] 
are multiplied according to (pb|). Inserting this solution into (|^) we obtain a new, equivalent representation for the 
conductivity 

Rea„^ = |^^5^«„(k)«Mk')E(-'^^)G-WG.(k){l- [Af,G,G.].}"'(k,k;k'-k) . (9) 

k,k' CTT 

For not too strong disorder the norm of the operator 11 A^'^.G^G,- 11 ^ 1. Hence the conductivity remains non-negative. 
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Note that only the non-local part with odd parity with respect to reflections in k and k' of the vertex 
A'^'' (k, k; k' — k) contributes to the conductivity. We obtain its leading asymptotic term if we use the represen- 
tation r = {1 - [A^'^GG] A^'' and solve it for A-^'', i.e., A'^'' = r{« [GGT] + Using the vertex T from (|) 
one finds in the order 0{l/d) 

A^''(ki, zi, k2, ^2; q) = A(zi, Z2) + (1 - A(zi, Z2)G(zi)G(z2))^ [r(ki, zi, k2, Z2; q) - r^''(k2 - ki; zi, Z2)] . (10) 

The irreducible vertex ( p^ together with the multiplication scheme (^) used in (|^) leads to an integral-equation 
representation of the conductivity. In the limit d — )■ 00 the momentum convolutions decouple. This fact helps us 
further simplify the expression for the conductivity. To derive the leading asymptotic contribution from the nonlocal 
part of A'^^ to the conductivity we have to calculate the momentum convolutions on the level of order 0{l/d) so that 
the velocities appear in squares and the momentum integrals do not vanish. 

In the following we resort to a hypercubic lattice where only the diagonal (longitudinal) conductivity remains. 
Keeping only the leading-order terms in the expansion of the denominator in (j^ we end up with a mean-field-like 
expression for the dc-conductivity 



where (i.jG,G,) = Af-> J^^vM'G,{k)Gr<k) and 



k,k' 



In the asymptotic limit d ^ 00 the irreducible vertex A'^'^ is determined from (|^) and (|T^) with a simplified momentum 

dependence of A*^'' via the quantity A'(k) — 2 J2 cos fc^, .111 We note that in the asymptotic limit d ^ 00 one has 

i/=i 

(A'^^) - 0(1) and (vlGaGr) 0{l/d). For {A"^^) = (as, for example, in the case of a k-independent Af^) the 
resulting expression for the conductivity reduces to the CPA-result 
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This is precisely the mean-fie ld co nductivity defined from the d ^ 00 limit, with Re crZ^ ^ 0{l/d) due to 



0(l/d). The denominator in (11a) then contains the leading asymptotic contribution from the vertex corrections to 
the conductivity. _ 

We note that the self-energy E and the vertex function A*^'' are connected via a Ward identity. VelickyE£l showed 
that the CPA self-energy E, (^), and the CPA-vertex A, (|^), fulfill the Ward-identity exactly. In our case, where we 
use the local CPA self-energy and the non-local part of the vertex function A*^'', (|l0|), or rather (A'J^^) from ( |llb| ), the 
Ward identity is fulfilled only asymptotically in the leading order of 1/d for both quantities. This is fully consistent 
with the spirit of the simplification we made in deriving the mean-field expression for the conductivity with vertex 
corrections (p]). 

Generally, in (^Tj) we have to perform integrals in momentum space that cannot be reduced to integrals over the 
density-of-states as it would be typical for mean-field theories. A reduction to an expression with integrals over the 
density-of-states is possible only if we resort to the leading contribution to (A'^^) in 1/d. This further simphfication 
yields 



Uvrd j 5 1 + t;l(G„Gr)A„r (1 - A^.G.G.) h„r(Gl)(Gi) - i^AGl)^ - irAGl)^] ' ^^^^ 
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{G^Gr)A,r (1 - A,,G,G,) [7.r(G2)(G2) - 7..(G2)2 - 7,,(G2)2] 

where (GcrGr) is defined as in (|ll| ) and "/err — i{Ef + iaO^jEp + «tO+). Here we used the fact that the velocity 
can be factorized, and (u^) = t*'^/2d on a hypercubic lattice in the high-dimensional limit. We immediately see that 
expression ( |l3|) is exact to 0(l/c?) for the conductivity and to 0{l/d^) concerning the vertex corrections. For the 
conductivity itself to be exact to 0{l/d^) an additional contribution in the numerator due to self-energy corrections 
of 0(l/c?) would have to be included. This, however, would not affect the vertex corrections in the leading order. We 
note once more that results obtained from a direct expansion of the conductivity or the self-energy in powers of l/d 
would not be guaranteed to be physical, e.g., the conductivity may become negative and the self-energy non-Herglotz. 
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Vertex corrections to the conductivity are particularly important in low dimensions (d = 1,2), where they lead to 
Anderson localization even for arbitrarily weak disorder. This is due to the diffusion pole in the two-particle propagator 
with energies from different complex half-planes. Eq. (pi]), which was derived in the limit of high dimensions, is not 
expected to be applicable in low dimensions, or to describe Anderson localization. Nevertheless, it is apparent from 
the form of the vertex in eq. ( [lO| ) that A^'^ includes a diffusion pole (in V^^) and thereby al so a Cooper pole describing 
coherent back scattering. For (—err) — +1 the k-integrals over the bare Cooper poles in ( |llb| ) diverge in dimensions 
d < 2. Hence d = 2 naturally appears as a critical dimension in our high dimensional approximation where the 
mean-field description must break down. 

In dimensions d > 2 the effect of the diffusion pole in on the electrical conductivity is not as strong in general 
and depends on the strength and type of disorder, band filling and the lattice structure. Below we demonstrate this 
trend with the help of two different disorder distributions. We put the Fermi energy into the band center and choose 
the dispersion of a d = oo Bethe lattice, with next-neighbor hopping amplitude i* = 1 to simplify the relation between 
the self-energy and the local propagator. First we choose the simplest disorder model, i.e. a percolation-type disorder 
distribution with Vi —0,oo occuring with probabilities 1 — x,x, respectively. The conductivity from (|l^) can now be 
evaluated explicitly. The result is 
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(14a) 



The conductivity (14a) remains positive everywhere and is a monotonically decreasing function of the disorder strength 
X. Fig. shows that the vertex corrections decrease the one-particle (i.e. CPA-) conductivity 
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FIG. 1. Mean-field conductivity with vertex corrections a = 
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1 as a function of the concentration x for percolation disorder in d - 



Another and more important example is the standard binary alloy with two values of the random potential Vi 
with equal probability, where we obtain 
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Fig. H shows that vertex corrections cause the conductivity to slightly increase with respect to the CPA result 

Ro .cPA,bi„ ^ (^\ 2(1 -A^) 
^47rdy A2(2- A2)' 



(15b) 



at weak and moderate disorder strength. Around A = l/\/2, the two-particle s catte ring in the vertex corrections 
result in a decrease o f th e CPA conductivity. For stronger disorder conductivity ( |l5a| ) remains positive and smaller 
than the CPA result (151) up to the split-band limit A = 1, where it vanishes. Both conductivities ( p^ monotonically 
decrease as functions of the disorder strength. However, the largest relative suppression of the CPA conductivity due 
to vertex corrections occurs at A = l/\/2 where the denominator of the second term in (15a), i.e., the terms with 
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(TT = 1 in ([l3|), diverges and changes sign. Note that if the denominators in (15a) were expanded only to the leading 
order in l/d, the conductivity would turn negative and hence become unphysical around A = l/v^. This explicitly 
demonstrates the importance of calculating the conductivity via approximations to the irreducible vertex A*^'', as 
proposed in this paper, instead of evaluating the conductivity corrections to the single-bubble diagram directly from 
the full vertex F. 




FIG. 2. 
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The results for the two models presented above demonstrate that the effect of vertex corrections in high dimensions 
is non-universal. Additional scatterings from the random potential contained in the vertex corrections can decrease 
as well as increase the single-bubble conductivity. When treated inappropriately vertex corrections can even lead to 
unphysical, negative results. The actual impact of vertex corrections in d > 2 depends on band structure, band filling, 
and disorder distribution and strength. 

In conclusion, we derived a mean-field expression for the vertex corrections to the electrical rfc-conductivity which 
becomes exact in the asymptotic limit of high lattice dimensions. To warrant the conductivity to be non-negative it 
was expressed as a functional of the irreducible vertex function in the electron- hole channel A^'^. We calculated the 
leading high-dimensional asymptotics of the non-local part of A"^'* and thereby derived an expression in closed form for 
the dc-conductivity including vertex corrections. Although this mean-field approximation does not describe Anderson 
localization, it goes systematically beyond the CPA. The result, Eq. (pT]), can be used as a mean-field formula for 
calculating the effects of vertex corrections to the electrical conductivity in, e.g., three-dimensional alloys, and may 
serve as a starting point for improved approximation schemes beyond the CPA limit. To include the essentials of 
the physics of Anderson localization one needs to improve the mean-field approximation for the irreducible vertex 
presented in this paper. A minimal requirement for this is a self-consistent theory for the non-local part of the vertex 
A'^'' which may, for example, be obtained from a parquet approximation. 
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